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■ Abstract 

C*"} ■ The Vafa-Witten arguments on the realization of parity and flavour 

symmetries in the QCD vacuum do not apply to two lattice regularizations 
of QCD which are able to reproduce the chiral anomaly: Wilson fermions 
and Ginsparg- Wilson fermions. We show here how using the last regu- 
larization one can get, from first principles, that the more standard order 
parameters for these symmetries take a vanishing vacuum expectation 
value for any number of flavours, if quarks have a non- vanishing mass. 
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Since the very beginning of the formulation of QCD as the gauge theory which 
describes the strong interaction between elementary particles, the understand- 



(•f-) ■ ing of the realization of symmetries in the QCD vacuum has been an important 



issue. During the 80's, Vafa and Witten gave arguments against spontaneous 
■ breaking of parity pQ and vector-like global symmetries [2J in vector-like theories 

such as QCD, and in principle, this seemed to settle the matter; however, these 
arguments were not as useful as expected. Indeed some years after the publica- 
tion of pQ , several articles appeared [21 HI El E] , calling into question the validity 
of the paper. The fact that the issue is still open twenty years after the publica- 
tion of the first paper is indicative of the complexity of the subject. Regarding 
vector-like symmetries such as flavour or baryon number conservation, it must 
be remarked that the Vafa-Witten theorem [5] is not applicable neither to the 
Ginsparg- Wilson regularization nor to one of the most used fermionic regulariza- 
tions on the lattice for QCD, i.e., Wilson fermions. In the Ginsparg- Wilson case 
the theorem does not apply because even if the integration measure is positive 
definite, the other essential ingredient in the proof in [2J, the anticommutation 
of the Dirac operator with 75 , is not realized. For the case of Wilson fermions 
neither of the two assumptions in [2] , positivity of the integration measure and 
anticommutation of the Dirac operator with 75, are fulfilled, the first of the two 
failing for an odd number of flavours. Indeed there exists a region of the pa- 
rameters space where parity and flavour symmetries are spontaneously broken: 
the well known Aoki phase [5] , and even a more complex phase structure for 
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lattice QCD with Wilson fermions has been recently suggested [TU]- In the 
end, a theoretical proof of the realization of symmetries of QCD is still lacking. 

The standard wisdom is that Vafa and Witten theorems fail, when applied to 
Wilson fermions, due to the existence of exceptional configurations which have 
a non- vanishing weight in the Aoki phase. Outside this phase, and in particular 
in the physical region near the continuum limit, the exceptional configurations 
would be suppressed and then parity and flavour symmetries would be restored 
in the QCD vacuum. Following this wisdom, it would be very convenient to 
choose a 'small eigenvalue free' regularization of QCD. It happens that Ginsparg- 
Wilson fermions fulfill this requirement. 

In this paper we will show how the probability distribution function (p.d.f.) 
of all the standard fermion bilinear order parameters for parity and flavour 
symmetry, as well as the p.d.f. of the topological charge density, are trivial in 
lattice QCD with Ginsparg- Wilson fermions of non- vanishing mass. The essen- 
tial ingredient in what concerns the proof of parity realization in the vacuum 
will be hermiticity. The rest of this article is organized as follows. In section 2 
we recall briefly the main features of Ginsparg- Wilson fermions and analyze the 
one-flavour case. Section 3 generalizes the results of the previous section to an 
arbitrary number of flavours, and shows how, contrary to the Wilson fermions 
regularization, the two flavour model has not an Aoki phase in any region of 
the parameter space. Section 4 contains a short discussion on the expectations 
about the realization of these symmetries in the chiral limit. In this limit the 
proof developed in the previous sections does not apply, and therefore we will 
argue following the standard wisdom. Hence, the content of that section is 
rather speculative. Last section contains our conclusions. 

2 Ginsparg-Wilson fermions on the lattice 

Years ago Ginsparg and Wilson (G-W) 12^ suggested, in order to avoid the 
Nielsen and Ninomiya non go theorem [131 114) and to preserve chiral symmetry 
on the lattice, to require the following condition for the inverse Dirac operator 



where a is the lattice spacing and R is a local operator. Accordingly D should 
satisfy, instead of the standard anticommutation relation of the continuum for- 
mulation, the Ginsparg-Wilson relation 



Fifteen years after this proposal, Hasenfratz [15] and Neuberger [16l [17] 
found that the fixed point action for QCD and the overlap fermions satisfy 
respectively the Ginsparg-Wilson relation, the last with R=l/2. Furthermore, 
Hasenfratz, Laliena and Niedermayer [IS] realized that Ginsparg-Wilson fermions 
have nice chiral properties, allowing us to establish an exact index theorem on 
the lattice. Indeed if we define a local density of topological charge as 



75 D- 1 +D- 1 j 5 =2aRj 5 , 



(1) 



7 5 £> + D75 = 2aDRj 5 D. 



(2) 



q(x) — aRTr (75-D (x, x)) 



(3) 



one gets a topological charge 



Q = aRTr{ lb D) 



(4) 
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which is a topological invariant integer that approaches the continuum topolog- 
ical charge in the continuum limit [T31 1201 I2T1 |2"21 I23] . 

Finally by replacing the G-W Dirac operator D + m by 

/ am\ , . 

A + m=[l —JD + m (5) 

in order to define an unsubtracted proper order parameter [24) 

^(l-^)v, (6) 

it is easy to see [25] that the G-W fermionic action possesses an exact sym- 
metry which is anomalous for the flavour singlet transformations, but exact for 
the flavour non-singlet case; a property which allows us to introduce also a 9 
parameter in the G-W action, as in the continuum. 

In this paper we will work with G-W fermions that realize relation ([2J with 
R = 1/2 (overlap fermions are a particular case of that), and with the massive 
Dirac operator (J5J associated to the unsubtracted chiral order parameter ^ . 

We derive in appendix |2] some properties of the spectrum of A + m and 
related operators that will be useful in the following. 



3 The one flavour case 

Let us start by considering the one flavour case. We will make use of the for- 
malism developed in [51111). First, we analyze the standard bilinear O = iip"f5ifi- 
Our basic tool is the function P(q), which is the generating function for the 
moments of O: P(0) = 1, P'(0) = i < O >, • • ■ , P"(G) = {i) n < O n >, where 
the average is the lattice gauge theory average with the full action. We can 
compute P(q) as [Sj 

det(A + m+f 75 ) det(g+£) m 
W { det(A + m) ; K det(JJ) ; 1 ' 

where V denotes the size of the matrix A (that is, up to a constant, the number 
of points in the lattice), H = 75 (A + to), and the averages are now taken with 
the effective gauge theory measure obtained after integrating out the fermion 
fields 

[(L4] e - SyM det(A + TO) (8) 
We denote by \ij the eigenvalues of H. Then we have for P(q) 

In appendix [X] we show that m < \fij\, so everything is well defined here as long 
as to > 0. Expanding the product we rewrite P(q) as 
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The sum in (|10[) is taken over all different combinations of k indices, so we can 
rearrange it as 



But we show in appendix [21 that the \i come in pairs ±/i, except maybe for the 
ones corresponding to chiral modes. Therefore most of the terms in ([TTj) cancel, 
and we are left with the contribution coming from the chiral modes: 

— >— = —— (n+ - n ~)+- (n' + - n'~) = ( 1 — 12 

V 3 H V \m y ' 2 K >) V2 J mV ' 

Putting everything together, we have, for non-zero mass 
v /, . ,„ , /,,\ k 



™-££(™-0(g)>+<4> 



(13) 



We see that in the thermodynamic limit the only dependence of P(q) on the 
gauge field is through the topological charge. From (TT3")) the moments of the 
probability distribution of O = iip^tp are related to the ones of the density of 
topological charge % through 

<(*»♦)-> -<-o-(= ! -i)"s^+<4) (»> 

All the odd moments vanish by symmetry. The first non-trivial moment is 
the second one 

( (l ^)V-(^-l) 2 i((|) 2 ) + 0(^) (15) 

At this point we require that vt/yy^tp be an hermitian operator. This is a sim- 
ple requirement, but with important consequences: The expectation value of 
the square of an hermitian operator must be positive, but from (I15p . in the 
thermodynamic limit it is manifestly negative. The only way to fulfill both 
requirements at the same time is for the second moment to vanish, 

Q N 2 



V^J{V)^ (16) 

But then the probability distribution of the density of topological charge, Q/V, 
must go to a delta at the origin 

JS.' 

and all the higher moments of both i^^ip and y vanish as well, therefore parity 
is not broken in lattice QCD with one flavour of Ginsparg- Wilson fermions, at 
least for the standard order parameter iipj^ip. 
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Let us consider now the case of the unsubstracted order parameter, 



#75 (l ~ ^ (18) 

We have 

de t (A + m +» 7 »(l-^)) 
det(A + m) 

Let's prove that only the zero modes of D contribute to (flT)|) . The matrix 
corresponding to the numerator is block-diagonal, and the contribution to the 
ratio coming from a pair of complex eigenvalues of D is of the form (from (|63[) 
in appendix lAl 

l-a£ (20) 

with \a\ < mr" 1 (l64l) . Therefore the contribution to P(g) corresponding to 
complex eigenvalues comes from the factor J\j fl — Q-7 yrj , where the product 
extends over all pairs of complex eigenvalues. Expanding this product, we have 
for the coefficient corresponding to q 2k 



yih a ^ ' ' ' a i* 

(ju—M 




m 2k 



(21) 



Thus the contribution from the complex eigenvalues is of order 1 + 0(y), that 
is, just 1 in the thermodynamic limit. 

The chiral modes of D with A = - contribute also 1 (lo r 5|) . The zero modes 
of on the other hand, give a non-trivial contribution: 



2 \ n+ 



If n+ <n~ (Q> 0), 



^n 1 -^ ^"^J (23) 

A similar expression is valid for n + > n~. We can argue essentially as in (12"TT) to 
see that the first factor in (123)) goes to 1 in the thermodynamic limit. Therefore 
we obtain the final result (valid for arbitrary values of n + , when — > cxd 



P(q) = (l- slgn (Q)^-) lQl (24) 
ml/ / 



All odd moments vanish as before because of (finite-volume) parity symmetry. 
The even moments also vanish because they are trivially related to the ones in 
(TT3l) , as can be seen easily by expanding ([24]) . In consequence we see that parity 
is not broken for the unsubstracted order parameter either. 

It is interesting to give a more physical argument that uses only the vanishing 
of the second moment. In fact if parity were spontaneously broken we would 
expect two degenerate vacua a and /3, since parity is a Zi symmetry. Let z a be 
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a complex number which give us the mean value of the pseudoscalar P in the 
a state 

(P)a = *»> 

then we have 

(P)p = ~z a - 

Since P 2 is parity invariant, it will take the same mean value in the two states 
and making use of the cluster property in each one of the two states we get 

(P 2 ) = \(P 2 ) a + \(P 2 ) l3 = *l 

but since (P 2 ) = 0, z a = 0. 

In conclusion we have shown rigorously, assuming hermiticity of i^^tf) and 
using standard properties of Ginsparg- Wilson fermions, that parity is not spon- 
taneously broken in the one-flavour model, at least for the standard order pa- 
rameters. 



4 The Np flavours case 

We study now the case of Np flavours, in general with different, non-zero masses. 
Most of the results from the previous section apply here as well with small 
modifications. The fermionic action is 

N F 

Y^tp a (A + m a )il} a (25) 

a=l 

The complete spectrum of the Dirac operator consists of Np copies of the single 
flavour spectrum, each of them calculated with the mass of the corresponding 
flavour. 

Let us consider the usual pseudoscalar order parameter for a single flavour /3, 
iippjsipp. The corresponding generating function P(q) can be computed easily: 

dct(A + m^f75) 
P{q) = { det(A + r4 )n * (26) 

The average is taken over the effective gauge theory with Np flavours, but only 
the (3 flavour appears within the averagqj- The calculation is identical to the 
one for the single flavour case and gives 



P(q) = (II ( 1 

3 V 



I) (27) 



The superindex on the eigenvalues indicate flavour, that is, belongs to the 
spectrum of 75 (A + mp). Equation (|2T|) is essentially the same as the one for 
one flavour, and the moments of P{q) still are, up to constants, the same as 
the moments of the density of topological charge Q/V. The only difference 
is that now the average is taken in the theory with Np flavours. This doesn't 
change any of the conclusions: all odd moments vanish by symmetry, the second 



1 Our fermionic determinant (and hence, its eigenvalues) will always refer to a single flavour. 
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moment vanish in the thermodynamic limit because of the hermiticity of i^^ij), 
the density of topological charge goes to a delta centered on the origin, and 
therefore all the higher moments vanish as well (when V — > oo). This is valid 
for any flavour separately, and so will be valid as well for any linear combination 
A a iip a j5ip a . The extension to the unsubstracted order parameter is trivial. 

Let us consider now the degenerate case, all flavours with equal nonzero 
masses. The action now enjoys flavour symmetry. As stated in the introduction 
of this paper, the Vafa-Witten theorem [5] for vector-like symmetries does not 
apply to G-W fermions because even if the integration measure is positive defi- 
nite, the Dirac operator does not anticommute with 75. However we can study 
this symmetry with the p.d.f. method as we did for parity. 

Consider first the case of two degenerate flavours and the standard order 
parameters ipT3ip and iip^Tsip. Proceeding as before we find for the first order 
parameter: 

p , v / det(A + m + »f)det(A + m- ? :f) x ,-pr ( q 2 V 

P t<l) = ( de t(A + m)2 )n * = ^{ 1 + W) Uf { ) 

where Xj are the eigenvalues of A + m. By the same argument we have use 
repeatedly before, P(q) — > 1 in the thermodynamic lirnifl 
Similarly for ii^jsTsij} we obtain 

p{q) = ( — — )n * = <ii )Nf m 

The same argument applies, and we have also that P(q) — > 1 in the thermody- 
namic limit. 

The calculations can be repeated easily for the unsubstracted operators, and 
the result is the same. We can then conclude that here is no Aoki phase for 
two flavours of Ginsparg- Wilson fermions with non-zero mass. This conclusion 
should not be surprising at all, for the spectrum of the Ginsparg- Wilson operator 
with a mass term is depleted of small eigenvalues (« ^). 

The extension of this proof to Np flavours is trivial, because it is easy to see 
that the preceding results apply to any pair of flavours, and therefore general 
flavour symmetry is realized. 



5 The chiral limit 

The results of the previous sections can not be extended in a straightforward 
way to QCD in the chiral limit. We lose the non-trivial lower bound on the 
spectrum of A + m, < m < |A|, and the Dirac operator has exact zero modes 
corresponding to gauge fields with non-trivial topology. 

As we can not get a definite conclusion on the realization of parity in the 
chiral limit from first principles, we will follow in this case the standard wis- 
dom. Let us consider QCD with one massless flavour. In this limit the action 
of the model has a new symmetry, the chiral U(l) symmetry, which is anoma- 
lous because the integration measure is not invariant under chiral U(l) global 

2 This conclusion rests only on the bound on the eigenvalues, and not on any other specific 
property of the Dirac operator. 
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transformations. The Jacobian associated to this change of variables introduce 
an extra term to the pure gauge action proportional to the topological charge 
of the gauge configuration, the 6- vacuum term, which allows one to understand 
the absence of a Goldstone boson in the model, but this fact generates the well 
known strong CP problem. All these features of QCD in the continuum formu- 
lation are well reproduced in lattice QCD with Ginsparg- Wilson fermions, as 
discussed at the beginning of this paper. 

Let's consider first the unsubstracted scalar order parameter, i/> (l — tp. 
The corresponding generating function in the chiral limit and for a generic value 
of 9 can be written as 

P(q) = 1+ — _i v_V U± (30) 

Zj 

with Z = J [dA] e ~ s YM e -*8Q det (D) and Q = n~ -n+. The contribution to the 
determinant in the numerator of (|30p coming from pairs of complex eigenvalues 
of D is easily computed and gives a factor 

M ? >=n(W-£(i-^)) (3D 

j 

where the product is taken over all different pairs of complex eigenvalues. Each 
chiral mode corresponding to an eigenvalue - contributes a factor of ^ to the de- 
terminant, and each chiral mode corresponding to a zero eigenvalue contributes 
a factor of ify. The normalization factor Z is of course calculated from the same 
expressions by setting q — 0. Therefore we can write P(q) as 

[dA]e-*"e-*«^J (JJ f (q) (32) 

The computation for the pseudoscalar order parameter iipj5 (l — ^p) i> is 
similar and the final result is 

f /l\ n ' ++n '~ + / ft \n++n- 

P(q)=Z- 1 J[dA]e- s ^e-^[^j C" 1 ) (f) M (33) 



n \ n + +n / ■ \ n +n 
2 \ / iq 



With the definition 

fs(l)=[^) [y) /o(?) (34) 

and denoting by Ps(q) and Pp(q) the generating functions for the scalar and 
pseudoscalar respectively, we can rewrite the above results in the following way 

P s (q) = Zg 1 J{dA}e- s ™e-^f s (q) (35) 

P P (q) = Zp 1 J [dA]e- s ^e~^f s {q){^ (36) 

From these expressions we see that only the Q = ±1 sectors contribute to 
the chiral condensate d (S), and that the dependence with the parameter 9 is 



3 In fact only a subset of the Q = ±1 configurations give a non-vanishing contribution, 
those with only one zero mode, that is, that verify n~ + n + = 1 as well as n~ — n + = ±1. 
Similar considerations apply to the other cases discussed in the text. 
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(S) oc cos(0). For the second moment (S 2 ), the only non- vanishing contributions 
come from the sectors with Q = and Q = ±2, and we can write 

(S 2 ) = A + A 2 (37) 

where Aq is the contribution coming from the Q = sector and A 2 the con- 
tribution coming from the Q = ±2 sector. Looking at the expression for the 
pseudoscalar condensate, all the odd moments vanish because of parity symme- 
try. For the second moment we find that 

(P 2 ) =A -A 2 (38) 

Looking at moments of higher order we would find an infinite set of relations. 

Since due to the chiral anomaly, strictly speaking, we do not have a new 
symmetry in the chiral limit of one flavour QCD, the standard wisdom is that 
the vacuum expectation value of the chiral order parameter will not vanish as 
m — > 0. Moreover, the absence of massless particles in the one flavour model 
suggests that the perturbation series in powers of m does not give rise to infrared 
divergences [27] . the free energy density is an ordinary Taylor series in m [261,127) ; 
and in what concerns the chiral condensate, the chiral and thermodynamical 
limits commute. 

On the other hand, the free energy density of the model atm^O and 8 = 
can be computed in the thermodynamical limit from the topologically trivial 
sector Q = • But since chiral symmetry in the Q = sector is not 

anomalous, it should be spontaneously broken in the topologically trivial sector 
if the chiral condensate takes a non vanishing value when approaching the chiral 
limit. In such a case, the value of the chiral condensate in the full theory and in 
the chiral limit will be related to the spectral density at the origin of eigenvalues 
of the Dirac operator of the topologically trivial sector by the well known Banks 
and Casher formula 28J 

(S) = -7rp(0) = S . 

This equation provide us with a non trivial relation between the value of the 
scalar condensate in the chiral limit in the full theory, which gets all its contri- 
bution from the Q — +1 sectors, and the spectral density of eigenvalues at the 
origin of the Dirac operator in the topologically trivial sector Q = 0. 

The scalar condensate is invariant under parity, and therefore in the full 
theory in the chiral limit, irrespective of the realization of parity in the vacuum 
we expect its probability distribution function to be a delta function <5(c — So). 
Therefore, we expect for the moments 

(S n ) = S» (39) 

For the second moment, this tells us that 

A + A 2 = Sg (40) 

But as previously stated, the standard wisdom tell us that the topologically 
trivial sector breaks spontaneously the chiral U(l) symmetry, and since A is 
exactly the second moment of -Ps(c) computed in the Q = sector, we should 
have 

A) = (41) 



9 



and therefore also 

A2 = (42) 

But this implies that the second moment of the pseudoscalar condensate van- 
ishes, 

(P 2 ) = (43) 

By a similar argument we can see that for higher (even) moments, (S 2n ) = Sq™ 
implies the vanishing of the corresponding pseudoscalar moment, (P 2n ) = 

Symmetry under parity is the only obvious reason for the vanishing of the 
pseudoscalar moments, and therefore the previous result strongly suggests that 
parity is also realized in QCD with one massless flavour. 

6 Conclusion 

Although the common lore on QCD symmetries states that parity and vector- 
like global symmetries remain unbroken, no sound theoretical proof of this hy- 
pothesis has ever been presented. The arguments given by Vafa and Witten 
against spontaneous breaking of these symmetries in [T| [5] were questioned by 
several groups [21 HI [SJ E] , and now there is agreement in the scientific commu- 
nity on the lack of a proof for parity realization in QCD. Concerning vector-like 
symmetries as flavour or baryon number conservation, it must be remarked that 
the staggered fermion discretization is the only known lattice regularization that 
fulfills the initial conditions of the Vafa- Witten theorem. Indeed the theorem 
is not applicable neither to the Ginsparg- Wilson regularization nor to Wilson's 
one, widely used for lattice QCD simulations. In the first case the theorem 
does not apply because even if the integration measure is positive definite, the 
Dirac operator does not anticommute with 75. In the second case neither of 
the two assumptions, positivity of the integration measure and anticommuta- 
tion of the Dirac operator with 75 are fulfilled. Indeed, there exists for Wilson 
fermions a region of parameter space where parity and flavour symmetries are 
spontaneously broken, the well known Aoki phase. 

On the other hand the p.d.f. formalism can be used to cast some light 
on the old aim of understanding the realization of symmetries of QCD from 
first principles. In fact, we have shown in this paper that some interesting 
conclusions appear when we apply the p.d.f. formalism to the Ginsparg- Wilson 
regularization. There, we see how the more standard order parameters for parity 
and flavour symmetries take a vanishing vacuum expectation value for a non- 
zero fermion mass. This is a major result that overcomes the difficulties found 

by mm- 
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A Spectrum of D and Related Operators 



We start with the Ginsparg- Wilson relation with R — 1/2, 

j 5 D + D l5 = aD l5 D. (44) 
We can also choose D such that 

l5 D l5 = (45) 

From (gU and ([35]) we obtain 

D + D^ = aDD^ = aD^D (46) 

Therefore D is a normal operator, and as such has a basis of orthonormal eigen- 
vectors. Also eigenvectors corresponding to different eigenvalues are necessarily 
orthogonal. From (|46p it is immediate to check that the operator V = 1 — aD is 
unitary, V'V = I. Therefore the spectrum of V lies in the unit circle with center 
in the origin, and the spectrum of D must then lie in the shifted and rescaled 
circle of radius - centered in the real axis at (-,0). The possible eigenvalues of 
D are therefore of the form 

A = -(l-e M ),a6R (47) 
a v ' 

We also have the identity 

A + A* = aAA* (48) 

Let v be an eigenvector of D with eigenvalue A, Dv = Av. Taking into 
account (l44l) 



-D75V = — 75-Dv + aDj^Dv = —A (75V + aDj^v) (49) 
Therefore using (|4T|) 

#(75v) = --^-( 7 5v) = A*( 75 v). (50) 
1 — aA 

Thus, if v is an eigenvector of D with eigenvalue A, then 75V is another eigen- 
vector with eigenvalue A* , and if A is not real then those two eigenvectors corre- 
spond to different eigenvalues and must be orthogonal. On the other hand if we 
restrict to the subspace corresponding to real eigenvalues, A = or A = -, 75 
and D commute, and therefore we can find a common basis of eigenvectors; in 
other words, we can find an orthonormal basis for which the eigenvectors of D 
corresponding to real eigenvalues are chiral. If we denote by n + (n~) the num- 
ber of eigenvectors of positive (negative) chirality in the subspace corresponding 
to A = 0, and similarly n' + (n'~ ) for the subspace corresponding to A = -, then 
Tr( 75 ) = and Q = f Tr( l5 D) imply 

n+ - rT = n'~ - n'+ (51) 
Q = rT -n+ (52) 

We denote by V the size of the matrix D. Then the density of topological 
charge, defined as is bounded in absolute value by 1, 



< 1. 
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The operator in the fermion action is 



A + m = (l - — ) D + m (53) 

Its spectrum is trivially related to the spectrum of D; if A are as before the 
eigenvalues of D, then the eigenvalues of (|53l) are (l — ^JA + m. They still 
lie in a circle with the center in the real axis, and the possible real eigenvalues 
are now m and — . We will always require that < m < -, then the operator 
([5"3")l preserves the position of the higher real eigenvalue [55] . 

We will also need the spectrum of H = 75 (A + m). It is easy to see that H 
is an hermitian operator, — H, and therefore has real spectrum We can 
calculate this spectrum by noting that the matrix 75 (A + m) is block diagonal 
in the basis of eigenvectors of D. If we denote by v A such and eigenvector with 
non-real eigenvalue A, we have 

ffv A = 75 (l ~ ^) Av A + mv A . =(m + \(l- ^)) v A » (54) 

ffv A , = 75 (l - ™) A*v A , + mv A = (m + A* (l - ^)) v A (55) 

We therefore have a 2 x 2 block 



m + A(l-2f) 

A* (1 - s*p) 



(56) 



The diagonalization of this block yields a pair of real eigenvalues ±/z with 

^ = m 2 + XX* (l - + m (A + A*) (l - ^) (57) 

For A real, let v A be an eigenvector of D of chirality x, that is, 75 v A = %v A . 
Then 

Hv x = 75 (l - ^) Av A + m X v x = (m + \(l- ^) ) *v A (58) 



That is, 



/' 



m + A(l-^)) X (59) 



More explicitly, we have /i = m and fj, = —m with degeneracy n + and n~ 
respectively, and similarly \i = #, ji = — # with degeneracy rt' + , n' — . 
Note that det(7s(A + to) = det(A + m), and therefore 

n*i=n^ ( 6 °) 

i 3 

From the above calculation we also obtain immediately a bound for fi at finite 
mass (this was remarked in |26j): 

pi 2 > m 2 (61) 

Let's consider now the operator 

A + to+--o< 5 ( 1- ^) (62) 



V V 2 
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We are interested in its determinant. Proceeding as before, we see that it is 
also block-diagonal in the basis of eigenvectors of D, and the contribution to 
the determinant coming from the block corresponding to a complex pair vx, 
is given by 



det 



1 



sp) A + m 



1 (-1 nA 
V \ 2 



m 



m 2 + 



am 



where we have used identity (|4"5|) . We also have the bound 



1 - £ 



AA* 



7/1 



1 



(am) 



AA* 



< 



A chiral mode A with chirality \ gives a contribution 
' am\ q ( aX\ 

1 --) x + m+ v( 1 -Y) x 



a 2 XX* 



(63) 



(64) 



(65) 



For A = - the contribution is iust -, whereas for a zero mode A = with 
chirality the contribution is m + Xy- 
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